Bright matter-wave soliton collisions at narrow barriers 
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We study fast-moving bright solitons in the focusing nonlinear Schrodinger equation perturbed by a narrow 
Gaussian potential barrier. In particular, we present a general and comprehensive analysis of the case where two 
fast-moving bright solitons collide at the location of the barrier. In the limiting case of a ^-function barrier, we 
use an analytic method to show that the relative norms of the outgoing waves depends sinusoidally on the relative 
phase of the incoming waves, and to determine whether one, or both, of the outgoing waves are bright solitons. 
We show using numerical simulations that this analytic result is valid in the high velocity limit: outside this 
limit nonlinear effects introduce a skew to the phase-dependence, which we quantify. Finally, we numerically 
explore the effects of introducing a finite-width Gaussian barrier. Our results are particularly relevant, as they 
can be used to describe a range of interferometry experiments using bright solitary matter-waves. 
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I. INTRODUCTION 

Bright solitary matter-waves are solitonlike dynamical 
excitations observed in atomic Bose-Einstein condensates 
(BECs) with attractive inter-atomic interactions [1-3]. They 
are solitonlike in the sense that they propagate without dis- 
persing [4], emerge largely unscathed from collisions with 
other bright solitary matter- waves and with external poten- 
tials [5, 6], and have center-of-mass trajectories which are 
well-described by effective particle models [7-9]. They de- 
rive these solitonlike properties from their analogousness 
to the bright soliton solutions of the focusing nonlinear 
Schrodinger equation (NLSE), to which the mean-field de- 
scription of an atomic BEC reduces in a homogeneous, quasi- 
one-dimensional (quasi- ID) limit. These bright soliton so- 
lutions of the ID focusing NLSE have been extensively ex- 
plored in nonlinear optics, both in the context of solitons in 
optical fibers [10-14] and as stable structures existing in ar- 
rays of coupled waveguides [ 15, 16] which are described by a 
discretized NLSE. Although the quasi- ID limit is experimen- 
tally challenging for attractive condensates [17], bright soli- 
tary matter-wave dynamics remain highly solitonlike outside 
this limit [3, 6]. Consequently, bright solitary matter-waves 
present an intriguing candidate system for future interferomet- 
ric devices [2, 6, 18-22]. 

A key component of a bright solitary matter-wave inter- 
ferometer is a mechanism to coherently split and recombine 
bright solitary matter-waves: the collision of a bright soli- 
tary wave with a narrow potential barrier is one way to create 
such a beamsplitter. Within a quasi- ID, mean-field descrip- 
tion of an atomic BEC, collisions of single solitary matter- 
waves with potential barriers and wells have been extensively 
studied [23-28], and sufficiently fast collisions with potential 
barriers have been shown to lead to the desired beamsplitting 
effect [27, 28]. When, in nonlinear optics, the soliton exists 
in an inhomogeneous array of discrete waveguides, the soli- 
ton can be reflected, split or captured at the position of the 
inhomogeneity [29-3 1 ] . This is equivalent, in the continuum 
limit of an infinite number of waveguides, to splitting a soli- 
ton in the GPE at a ^-function [29] — a phenomenon which 
has been called the "optical axe" [14]. Such splitting has 



been considered in the context of soliton molecule formation 
[21], within a mean-field description, and also in the context 
of many-body quantum mechanical descriptions: in the lat- 
ter it has been demonstrated that macroscopic quantum super- 
positions of solitary waves could be created, offering intrigu- 
ing possibilities for future atom interferometry experiments 
[19, 20]. Recently Martin and Ruostekoski, in Ref. [22], con- 
sidered an interferometer using a narrow potential barrier as 
a beamsplitter for harmonically trapped solitary waves, based 
on the particular configuration of a recent experiment [32]. In 
particular this work demonstrated that such a potential barrier 
can also be used to recombine solitary waves, by arranging 
for them to collide at the location of the barrier. In such col- 
lisions, the relative norms of the two outgoing solitary waves 
was shown to be governed by the phase difference A between 
the incoming ones. In the mean-field description the rela- 
tive norms of the outgoing waves exhibit enhanced sensitivity 
to small variations in the phase A; however, a simulation of 
the same system including quantum noise, via the truncated 
Wigner method [33], showed increased number fluctuations 
that ultimately negated this enhancement [22]. 

In this paper we consider the focusing NLSE perturbed 
by a narrow, Gaussian potential barrier of the form V(x) = 
qe~ x ~l 2(T ~ I y/2ncr, and investigate the dynamics of two fast- 
moving bright solitons which collide at the location of the 
barrier. We investigate such collisions for the general initial 
condition [Fig. 1(a)] 
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with b > 0. For large xq this approximates an exact two- 
soliton solution comprising two bright solitons with unequal 
norms, 1/(1 + b) and b/(l + b), oppositely directed and equal 
velocities, ±v, and relative phase A [12]. 

By examining such collisions for general b, A, q, and <x 
we give a detailed explanation of the nonlinear recombination 
which occurs after the solitons collide at the potential barrier 
at time t\ - 2xq/v, and are recombined into left- and right- 
travelling waves in a phase-sensitive way. This general and 
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FIG. 1 . (Color online) (a) Schematic of the collisions we consider: 
two bright solitons [one in the case b = (a)(i)] (solid lines) col- 
lide at a narrow Gaussian potential barrier (dashed line). The norms 
of the two outgoing waves are nonlinearly dependent on the relative 
phase A between the solitons, as illustrated in (b) for equal-amplitude 
solitons [the case b = 1 (a)(iii)]; solid red (dashed blue) lines indi- 
cate the outgoing wave in the negative (positive) x domain. Here the 
soliton velocity is v = 2 and the barrier width is characterised by 
cr = 0.28 



comprehensive treatment of two-soliton collisions at a bar- 
rier constitutes the main result of the paper. For the case of 
solitons of equal size (as reported in Ref. [22]) we illustrate 
this phase dependence in Fig. 1(b). In this paper we present 
an analytic description of the recombination for the general 
two-soliton case (b > 0) in the limit of a <5-function barrier 
(cr — > 0). This description is derived from an exact description 
of the single-soliton case (b = 0) in the same limit [27, 28]. 
We compare this to numerical simulations, and find the ana- 
lytic description is exact in the limit of high velocity. In ad- 
dition to yielding useful predictions for the relative norms of 
the recombined waves, this analytic method allows us to esti- 
mate whether one, or both, of the outgoing waves are bright 
solitons. We also numerically investigate the case of a Gaus- 
sian barrier, cr > 0. Particular cases of interest are b — 
(b — > oo) — corresponding to a single soliton — and b — 1 — 
corresponding to equal-sized solitons; these correspond, re- 
spectively, to the splitting and recombination stages of a bright 
solitary wave interferometer. While in the context of atomic 
BECs the NLSE represents a quasi- ID condensate with tight 
radial trapping and either zero or very weak axial trapping 
(e.g., a periodic "ring" trap [34], or a waveguide or weak har- 
monic trap [17]), we emphasise that the equation we study 
here remains general and could also be used to describe simi- 
lar systems in, e.g., nonlinear optics. However, as a particular 
example, our analysis directly allows us to understand the op- 
eration of a bright solitary wave interferometer in a ring trap, 
illustrated schematically in Fig. 5(a). 

The paper is structured as follows: in Section II we de- 
rive the NLSE, perturbed by a narrow Gaussian barrier, in the 
context of an attractively-interacting atomic BEC. The subse- 
quent sections comprise our analysis of the collisions given by 
initial condition Eq. (1). We consider first the single-soliton 
case {b = 0), for (5-function (Section III A) and Gaussian 
(Section III B) barriers, and subsequently the two-soliton case 
(b > 0), again for 5-function (Section IV A) and Gaussian 
(Section IV B) barriers. In section Section V we conclude by 
interpreting our results in the context of current and future 



atomic BEC experiments. 



II. PHYSICAL SYSTEM 

In general a weakly interacting atomic BEC, in the limit of 
zero temperature, can be described by the 3D Gross-Pi taevskii 
equation (GPE) [35]: 



ih 



d¥(r) 
dt 



2m 



V 2 + V t ™(r) + V„ t (r)+ ftD pi'(r) 



T(r). 



(2) 

Here gm - 4-nn 2 a s N/m, and N, m and a s are the atom num- 
ber, mass, and s-wave scattering length respectively. For at- 
tractive inter-atomic interactions a s < 0. The wave function, 
*P, is normalised to 1. The potential V t rap(r) = mu> 2 (y 2 + z 2 )/2 
represents the trapping potential, which we take to be a cylin- 
drically symmetric waveguide; such a configuration is approx- 
imately achieved in an atomic waveguide trap, or in a toroidal 
"ring" trap [34] which also introduces periodicity in x. 

By increasing the radial trapping one can reach a quasi- 
1D regime, as defined in detail in Ref. [17], where the ra- 
dial trapping is tight, but not such that the scattering is no 
longer 3D [a s <K (ft/m&v) 1 ' 2 ]. In this regime we can sep- 
arate the radial and axial dynamics with the ansatz *P(r) = 
^mixXmajr/ntiy^expi-mcJrly 2 + z 2 ]/2h). After factoring 
out global phases associated with the radial harmonic ground 
state energies, this yields the quasi- ID GPE 
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The nonlinearity is quantified by g — 2nu> r a s . We model the 
external potential as 
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This can be generated by an off-resonant Gaussian light sheet 
propagating in the z direction with 1 /e 2 radii x r and y r (y r » 
x r ). In this case A = u)l-u>o is the detuning of the light sheet's 
frequency u>l from the optical transition frequency cot), and Q 
is the Rabi frequency at the centre of the light sheet [36]. 

Working in "soliton units" — position units of ft 2 /mgN, 
time units of t? /mg 2 N 2 , and energy units of mg 2 N 2 /n 2 [17] 
— yields the dimensionless, quasi- ID GPE 
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where the dimensionless wave function is \p — WYyoJ yfmgN, 
the normalised barrier width is cr — (n 2 /2mgN)x r and the bar- 
rier strength is given by 



jc r Q 2 V2^r 
32io r a s NA' 
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III. ONE-SOLITON SPLITTING ON A NARROW 
BARRIER (b = 0) 

A. <5-function barrier (cr — » 0) 

In this section we examine the splitting of a single bright 
soliton {b = 0) on a <5-function barrier. The assumption of a 
<5-function barrier facilitates an analytic treatment and is valid 
for narrow barriers with cr — » 0. A detailed analytic treatment 
single-bright-soliton splitting on such a barrier is given by 
Holmer, Marzuola and Zworski in Ref. [27]. Here we briefly 
restate two key results of Ref. [27] within our own notation. 

Firstly, the transmission coefficient for a fast-moving bright 
soliton splitting on a 5-function barrier is approximately equal 
to the transmission coefficient for plane waves incident on 
an identical (5-function barrier in linear quantum mechanics, 
T q (v), given by 



T q (v) = \t q (v)\ 2 = 
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v 2 + q 2 



1 + a 2 



(7) 



Here, t q {v) is the transmission amplitude associated with a 5- 
function barrier in linear quantum mechanics, and the soliton 
velocity v plays a role analogous to the wavenumber of the in- 
cident wave. The transmission and reflection amplitudes t q (y) 
and r q (v), are defined as 
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IV 
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The quantity a characterises the transmission in the linear 
case, and hence the transmission of bright solitons in the high 
velocity limit. The exact relation between T q (v) and the actual 
transmission coefficient for the incident bright soliton, 



T s q { 
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\if/(x, t)\ 2 dx, 



(9) 



is determined in Ref. [27] to be 
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as v — > oo, 



(10) 



provided that the initial offset is xo < —v n and a - qjv 
is fixed. Here, 77 is a parameter linked to the duration for 
which the soliton interacts with the barrier, and must satisfy 
2/3 < T] < 1. The brevity of this duration for a fast-moving 
bright soliton, which allows one to treat the splitting as a lin- 
ear process, is fundamental to the proof of the above result 
[27]. The error term in Eq. (10) is minimized for brief colli- 
sions (77 — > 1), in which case it decays with velocity as v~ 112 . 
Secondly, it is also determined in Ref. [27] that the outgoing 
waves resulting from splitting a bright soliton on a 5-function 
barrier are composed of either one, or two, bright solitons, and 
a time-decaying radiation term. This is significant, as previ- 
ously the transmitted and reflected waves were considered to 
be only 'soliton-like' [26, 28]. The resulting bright solitons 
are described, for high velocity, by 

-\ X0 \/vT l ' 2 ) + O(v 1 -^ 2 ) 
(11) 



i//(x, t) - if/ T (x, t) + if/ R (x, t) + 0{[t 



where 

i/f T (x, t) = & 1 e i(xv+lAT - v2]l/2) A T sech(A T [x -x - tv]), 
ifr R (x, t) = e iVs e'' ( -' n ' +[AK - l ' 2l ' /2) A s sech(A s [jt: + x + tv]). 

The amplitudes of the transmitted and reflected solitons are 
given by 

A T = max(0, 2\t q (v)\ - 1) and A R = max(0, 2\r q (v)\ - 1); (12) 

in the case that At (A r ) is equal to zero, the transmitted (re- 
flected) outgoing wave does not contain a soliton, but only 
radiation. More generally, the inequalities Aj < T q (v) and 
A R < 1 - T q (v) hold. The phases imparted by the splitting 
process are defined by 



(fir = arg(f ? (v)) + ip Q {\t q {v)\) + [1 - A 2 ]|x |/2v, 
tp R = arg(r 9 (v)) + ip Q (\r q (v)\) + [1 - A 2 ,]|x |/2v, 



(13) 
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B. Gaussian barriers (cr > 0) 

We now analyse, numerically, the bright soliton splitting 
process at a Gaussian barrier. Our numerical simulations use a 
Fourier pseudospectral split-step method with a periodic grid. 
We ensure that grid size and spacing are chosen such that the 
bright solitons are well separated and the effects of the peri- 
odicity are negligible. 

Our initial condition takes the form 



^(x) = -e' VA sech([x-* ]/2), 



(15) 



where xq < 0. Figure 2(a) shows the transmission coefficient 
T q (v) obtained from numerical simulations of a single bright 
soliton splitting on a Gaussian barrier with width cr = 0.1, and 
with a = qjv - 0.6, 0.8, 1.0, 1.2, and 1.4. In our numerics 
we define T q (v) by the integral of </>(x, t{) over the positive x 
domain, 
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Here t\ - 2\xq\Iv, such that at this time an unimpeded bright 
soliton would have reached the point x - +\xq\; at this time the 
outgoing waves are well-separated. The results are compara- 
ble to the 5-function barrier case explored in Ref. [27] and the 
previous section. 

Figure 2(a) shows that as a increases so does the discrep- 
ancy between the asymptotic <5-function limit and T'(y). This 
can be understood by considering how the strength of the bar- 
rier compares to the (particle-like) kinetic energy of the soli- 
ton v 2 12. In the region where the strength of the barrier is 




considered, boils down to an argument that we must have 
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FIG. 2. (Color online) (a) Plot of numerically obtained bright soliton 
transmission, T°(v), as a function of velocity v for a range of fixed 
a = q/v and a narrow Gaussian barrier with width cr = 0.1. Dashed 
lines show the transmission through a (^-function in the linear regime 
for the same range of a. (b) Numerically obtained bright soliton 
transmission for a = 1 and with a range of barrier widths <x. 



greater than the soliton's kinetic energy the wave function de- 
cays, reducing transmission. By equating these two values, 
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to be definitely out of the classical transmission regime. From 
Eq. (19) it is apparent that for satisfactorily large v we will 
always enter the classical transmission regime for any given 
finite Gaussian barrier. This regime cannot be retrieved in the 
5-function case. 

The comparison to the 5-function case is valid in the quan- 
tum transmission regime, where the velocity is low enough 
(for a given q, cr) that the soliton cannot classically pass 
through the barrier and must tunnel through instead. For ex- 
ample, this is true when 0.5 < v < 2 and cr < 0.28 [Fig. 2.] 
Within the quantum transmission regime [Eq. (19)] the 6- 
function limit of 0.5 is reached (from below) by reducing cr. 
This allows for larger values of v, as is consistent with Holmer 
and Marzuola's work in [27] where results are general for any 
v > 1 (and so is in the high velocity regime). 

Figure 2(b) shows that the transmission approaches the an- 
alytic prediction for a 5-function barrier as the barrier width 
cr tends to zero. This confirms that the analytic expressions 
given in Ref. [27] and the previous section for the 5-function 
barrier can be quantitatively useful for realistic Gaussian bar- 
rier widths. For example, Fig. 2(b) indicates the analytic pre- 
diction is reasonably quantitatively accurate for cr < 0.28 in 
soliton units. For a condensate of 85 Rb and using typical ex- 
perimental parameters of N ~ 6 x 10 3 atoms, a s ~ 5«o (the 
Bohr radius) and a> r ~ 17Hz this translates to a splitting beam 
with a full width at half maximum of ~ 9 /mi. These param- 
eters are consistent with the experimental setup in [3]. For 
a similarly sized condensate of 7 Li atoms tuned to a similar 
scattering length this width becomes ~ 2 fim. This parameter 
regime is consistent with [1] appart from the radial trapping 
frequency, which we reduced from 27r x 7 10Hz to 2n x 200Hz. 



we determine that the distance over which the wave function 
decays, Xd, is described by 



IV. TWO-SOLITON COLLISIONS AT NARROW 
BARRIERS (b > 0) 



x 2 rl = 2cr 2 ln\J-^- 



2 a 



(18) 



It is clear that, for a given v and cr, as we increase a (by in- 
creasing q) we increase Xd- This is inconsistent with the as- 
sumption of a brief barrier-soli ton interaction period, which is 
required in the delta function case of soliton splitting. This in- 
sconsistency causes an increase in the attenuation of the wave 
function, reducing transmission. 

We show the computed dependence of the transmission on 
the barrier width cr in Fig. 2(b). These computations were 
carried out with a — 1 . For wider barriers or in the higher ve- 
locity range, where the peak height of the potential is less than 
the (particle-like) kinetic energy v 2 /2 of the incident soliton, 
the amount of transmission is greatly increased. This illus- 
trates the classical transmission regime where the soliton sim- 
ply passes through the potential, and, for the Gaussian barriers 



A. Analytic treatment for ^-function barrier (<x — > 0) 

We now give an approximate analytical description of the 
dynamics of two fast-moving bright solitons colliding at a 6- 
function barrier, which we subsequently compare to numer- 
ical simulations in order to give a fuller picture of the real 
dynamics that we might expect to see in an experiment. This 
analysis stops short of the full analytic rigor used in [27] but 
is consistent within its assumptions of linearity. As previously 
stated, during the time over which one bright soliton inter- 
acts with the potential we can describe the system as linear 
[27]. Here we extend this argument to a scenario in which two 
bright solitons collide at a 5-function potential, as described 
by the equation 
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'2dx 2 



q6(x)-\if,(x,i)Y 
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(20) 



and initial condition 

2 + 2b \2 + 2b) 



(21) 



[Eq. (1)]. We achieve this by making use of the second re- 
sult of Ref. [27], which we apply to the positive and negative 
domain bright soli tons, if/ + and i^_, separately, before taking 
a linear combination of the results. This means that at some 
time \xo\/v < t < v~ n + \xo\/v after the barrier collision the 
solution can be written as a sum of four sech profiles, two in 
each of the positive and negative domains; 

i//(x, t) = ip +T (x, t) + t//- R (x, t) + ifr_ T (x, t) + il/ +R (x, t). (22) 

Here 1/1+7- denotes the bright soliton transmitted to the negative 
domain which originated in the positive domain, t//_s denotes 
the bright soliton originating from and reflected back into the 
negative domain, and so on. In this scheme 

tf/ +T (x, t) = e'' w+r+ ^ +A) A +7 -sech (A +T [x - x Q + tv]) , 
ij, +R (x, t) = e ;( ^ + ^ +A) A +s sech(A +R [x + x - tv]) , 
il/- T (x, t) = e'^+^A-rsech (A_ r [x + x - tv]) , 
ifr. R (x, t) = e'^^U-ssech (A_« [x - x + tv]) . 

(23) 

Two phase factors appear above; the (/> ±R /t are those associ- 
ated with the standard soliton solution and are given by 



b ±T = +vx + [A 2 ±T - v 2 ] t/2, 
b ±R = ±vx + [A 2 ±R - v 2 ] t/2. 



(24) 



The <p ±R /T factors are imparted by the collision, and are de- 
scribed by 

<P±t = [l -A^] \x \/(+2v) + arg(f,(v)) + <p Qt q (+v)\), 

<P±r = [l -A 2 ±R ] \x \/(+2v) + arg(r 9 (v)) + <fio{\r q (+v)\). 

(25) 

With b = 1, barrier height q - v, and fast-moving solitons 
(v large) both initial bright solitons are split equally, such that 
the amplitudes A ±R /t are all equal and global phases can be 
dropped. In this case Eq. (22) simplifies dramatically, and 
shortly after the collision can be written as 



i//(x, t) = i//++(x, t) + if/—(x, t)), 
\j/ ++ (x,t) = P + (A)f + (x,t), 
t-(x,t) = P.(A)f.(x,i), 



(26) 



where the terms 



p (/\\ - _ | e 'arg(f,(?)) + g'[arg(f,(?))+A] 1 
p tfc\ - _ / e 'arg('?(?)) + gi[ais(r,(?))+A] J 



(27) 



contain information about the constructive and destructive in- 
terference between the transmitted and reflected waves. It 
should be noted that this treatment allows us to infer the bright 
soliton interactions, but does not give us a complete solution; 
the terms /+ and /I contain information about the outgoing 
wave profiles. By taking a linear superposition of the resul- 
tant bright solitons we initially obtain a sech profile which is 
not a single-soliton solution. However, in subsequent nonlin- 
ear evolution this profile returns to a soliton profile to within 
a known error, as documented in Appendix B of Ref. [27]. 

At a suitably large time after the collision, when the soli- 
tons have again separated to the extent that they are again ef- 
fectively independent, inspection of |i/<i 2 shows that the bright 
solitons are modulated by the factors 



|P_(A)| 2 =i[l-sin(A)], 
|P + (A)| 2 = i[l+sin(A)]. 



(28) 



These factors determine the norm of the outgoing waves in the 
positive and negative domains, defined by 



Jr»±oo 
\i[f(x,t)\ 2 dx = \P ± (A)\ 2 . 
o 



(29) 



Within the analytic approach presented here T ± are functions 
of A alone. It should be noted that the symmetry of the initial 
condition and linear interaction means that the phase interac- 
tions apply to both the transmitted and reflected bright solitons 
and the radiation terms. As a result the quantity T ± scribes 
the total density in the positive and negative domains, not just 
the respective bright solitons. For suitably high incident ve- 
locities this radiation becomes negligible, in accordance with 
Eq.(ll). 



B. Numerical treatment for (5-function and Gaussian barriers 
(general <r) 

In Fig. 3 we present results of numerical simulations of fast 
(v > 1) bright soliton collisions at both 5-function 1 and Gaus- 
sian barriers. The norms of the outgoing waves, defined in our 
numerics by 
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\ij/(x, t\)\~dx, 



(30) 



agree qualitatively with the predictions of our analytic treat- 
ment, but with a noticeable skew in the predicted sinusoid. 
This skew is also visible in the results for the Gaussian bar- 
rier case shown in Fig. 1. We parametrise this skew by e and 
describe the norms of the outgoing waves, T ± , as 



7V = 



1 ± sin(A + e) 



(31) 



1 Within our Fourier pseudospectral method a 5-function barrier can be im- 
plemented with high accuracy in momentum space using the approach out- 
lined in Ref. [37]. 
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FIG. 3. (Color online) Phase skew of numerical results with respect 
to analytic prediction for equal-size (b = 1) bright soliton collisons 
at a narrow barrier, (a) Numerically obtained data showing the de- 
pendence of the norm of the outgoing wave in the positive domain, 
T + , after Gaussian barrier collision, on phase difference A. Shown 
here are interference curves for solitons moving with velocity v = 1 
(fuchsia, solid), 2 (light blue, long dashed), 3 (red, short dashed), 4 
(green, dotted) and 5 (dark blue, dot-dashed). The width of the bar- 
rier was <x = 0.14. (b) (see (a)) Collisions at a ^-function barrier. 
Notice the qualitatively identical form of the curves, illustrating that 
both (5-function and Gaussian barriers exhibit the same skew, and so 
both undergo the same non-linear effects, (c) Numerically obtained 
data showing the phase perturbation e [Eq. (31)] due to non-linear 
effects in a soliton collision at a (5-function barrier. Shown here, in 
order of descending amplitude, are the skewness parameters (e) of 
the interference curves for solitons moving with velocity v = 1,2, 
3, 4, and 5 colliding at a <5-function barrier, (a)-(c) are the upper 
left quarters of the full data set; the plots are both symmetric about 
the lower and right hand axes. All results shown are calculated for 
a = q/v = 1. 



This skewness parameter is less pronounced for increasing ve- 
locities, i.e., 



lim max(e) = 0. 



(32) 



The presence of the skew in simulations with both Gaus- 
sian and ^-function barriers rules out any explanation in terms 
of the barrier structure. However, it is well known that when 
solitons collide in the absence of a barrier they induce a small 



phase and position shift in one another [12, 21, 38]. We pro- 
pose that the skew is a result of interactions between the soli- 
tons while approaching the barrier; more fundamentally, this 
is a result of the condition of a brief interaction not being fully 
satisfied. For instance, from initial condition Eq. (1) the phase 
(ip' t ) and position (xj) shift on the left hand soliton are given 
by 



2x\ 
1 + b 



+ iip'/ - 2 In 



v + ; 



v + i[(l - b)/(l + b)] 



(33) 



In the case of equal amplitudes and velocities total phase dif- 
ference reduces to </>' = ±4 arctan( 1 /v) or, in the limit of high 
velocity, tp' « +4/v. In our scenario only part of this phase- 
shift can occur before the solitons enter the linear regime, and 
so we expect that our skewness parameter e will be some frac- 
tion of if' . What we have observed from our numerics is that 
e oscillates with A but the maximum value is e max » </?'/8. 
This is consistent with the behaviour we observe in the high 
velocity limit. 

It should also be noted that the interference effect is present 
in collisions between solitons of differing amplitudes. By tak- 
ing b - e 8 we see that there is still interference between the 
transmitted positive and reflected negative bright solitons (and 
vice versa) [Fig. 4]. Along the line/? = 0, where the ampli- 
tudes of the incoming bright solitons are equal, we can clearly 
see a sinusoidal dependence on A. For nonzero f3 there is still a 
notable dependence on the incoming phase difference, but this 
effect is soon washed out if the difference in initial amplitudes 
becomes too large. It is true, however, that the solitons do 
not have to be of similar size to constructively or destructively 
interfere. 

The black (white) contour on Fig. 4 shows where the fi- 
nal population in the positive (negative) domain was not great 
enough, after interference, for the final aggregation to form a 
soliton. This is determined by treating the total positive (neg- 
ative) domain population as the transmission (reflection) coef- 
ficient in equation (12). From Eq. (12) we see that \t g \ and \r g \ 
must both be > 0.5 to get two outgoing solitons. We determine 
\t t] \ and \r q \ numerically as 



V^ 



if: 



\ip\ 2 dx. 



(34) 



As such, the white contour marks where T - 0.25 (T + - 
0.75) and the black contour marks where T + = 0.25 (T_ = 
0.75). 



V. APPLICATIONS AND CONCLUSIONS 

As stated in Section I, an important aspect of our analysis 
is that is directly allows us to describe the operation of the 
bright solitary matter- wave Mach-Zender interferometer in a 
ring trap shown in Fig. 1(c). In the quasi-lD limit we have 
considered in this paper, such an interferometer can be de- 
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FIG. 4. (Color online) Numerically computed transmission coeffi- 
cient T + illustrating the interference between solitons of different ini- 
tial amplitudes (b = eP) colliding at a (5-function barrier. Even in the 
case of a large difference in initial amplitude (large |/3|) there is still 
interference between the solitons. The contour lines show the bound- 
ary between having one (interior regions) and two (exterior region) 
outgoing bright solitons in the analytic treatment [see Eq. (12)]. All 
results shown are calculated for a = q/v = 1 and v = 5. 



cated on opposite sides of the ring trap. Our analysis can 
be applied to understanding such an interferometer by split- 
ting its operation into the following three stages (illustrated in 
Fig. 5): 

Firstly, a single initial bright soliton is incident on the bar- 
rier at x - (Fig. 5(a)(i)), at high velocity. Assuming an 
initial displacement of x$ - L/4 is sufficient for the soliton to 
be well-separated from the barrier, and a barrier height q - v 
(such that a - 1), the analysis of Section III applies, with 
b — 0. Hence, we obtain two equal-sized outgoing bright soli- 
tons, with relative phase tt/2. 

Secondly, these bright solitons propagate without disper- 
sion in opposite directions around the ring. We assume that 
the soliton in the positive x-domain picks up an additional 
phase shift A due to the effects of whatever interaction the 
interferometer is measuring [Fig. 5(a)(ii)]. 

Thirdly, these bright solitons collide at the barrier at x = 
L/2. Here, the analysis of Section IV applies, with b = 1, 
a = 1, and A — > A - n/2 [Fig. 5(a)(iii)]. In our analytic 
treatment, this means that the norms of the outgoing waves 



1 ± cos(A) 



(36) 
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FIG. 5. (Color online) (a) Schematic of a ring-trap interferometer; 
(i) an incoming bright soliton (dotted, green) is split into two equal- 
amplitude solitons at the first narrow barrier (dashed, black) gain- 
ing relative phase ~ tt/2; (ii) these solitons (solid, cyan) propagate 
around the ring, accumulating an additional relative phase difference 
A; (iii) at the second narrow barrier (dashed, black) these solitons 
are recombined into outgoing waves [dot-dashed, red (blue) for pos- 
itive (negative) x domain]. The norms of the two outgoing waves are 
shown as a function of A in (b), and illustrate the shift by ~ n/2 with 
respect to Fig. 1(b) caused by the initial splitting. Here the soliton 
velocity is v = 2 and the barrier width is characterised by a = 0.28 



1 d 



scribed by 
d\p(x) 



where x e (-L/2, L/2] is now a periodic coordinate. In 
Eq. (35) there are two narrow Gaussian potential barriers, lo- 



2 8x^ cr^ 

+ e- ( *- L/2W ]H«K*)l 2 W), (35) 



The predicted and computed A-dependence of T ± is shown in 
Fig. 5(b). The skew with respect to the analytic prediction we 
quantify in Section IV corresponds to the nonlinear enhance- 
ment of the phase-dependence reported in Ref. [22]. 

To conclude, we have presented a general and detailed anal- 
ysis of the collision of two fast-moving bright solitons at a 
narrow potential barrier in the NLSE. We have developed an 
analytic treatment of this problem, based on the assumption of 
a (5-function potential and short collision times. Our numerical 
simulations of the same problem reveal that this analytic treat- 
ment is quantitatively accurate in the limit of narrow barriers 
and fast solitons as described in Section IIIB. At realistic soli- 
ton speeds and barrier widths, however, our numerical results 
differ from the analytic prediction; we have quantified this in 
terms of the phase-skew e. Our analytic treatment also pro- 
vides an estimate of the regimes in which the outgoing waves 
contain solitons. One important application of our analysis is 
describing the operation of a bright solitary matter-wave inter- 
ferometer in a ring trap. However, we stress that our analysis 
remains general, and could potentially be used to describe a 
range of possible interferometry experiments, either in bright 
solitary matter-waves or other physical systems. 
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